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Abstract. In this paper we continue to explore the connection between tensor algebras 
and displacement structure. We focus on recursive orthonormalization and we develop an 
analogue of the Szego type theory of orthogonal polynomials in the unit circle for several 
noncommuting variables. Thus, we obtain the recurrence equations and Christoffel-Darboux 
formulas for Szego polynomials in several noncommuting variables, as well as a Favard type 
result. Also we continue to study a Szego type kernel for the A^-dimensional unit ball of an 
infinite dimensional Hilbert space. 
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1. Introduction 

In the first part of this paper, 0, we explored the connection between tensor algebras 
and displacement structure. The displacement structure theory was initiated in [[L3[] 
as a recursive factorization theory for matrices whose implicit structure is encoded by 
a so-called displacement equation. This has been useful in several directions including 
constrained and unconstrained rational interpolation, maximum entropy, inverse scat- 
tering, if °°-control, signal detection, digital filter design, nonlinear Riccati equations, 
certain Fredholm and Wiener- Hopf equations, etc., see | Iifl . Aspects of the Szego the- 



ory can be also revealed within the displacement structure theory. Our main goal is to 
develop an analogue for polynomials in several noncommuting variables of the Szego 
theory of orthogonal polynomials on the unit circle. An analogue of the Szego theory 
of orthogonal polynomials on the real line is being developed in the companion paper 

0- 

The paper is organized as follows. In Section 2 we review notation and several re- 
sults from |J. In this way, this paper can be read independently of 0. In Section 3 
we introduce orthogonal polynomials in several noncommuting variables associated to 
certain representations of the free semigroup and discuss their algebraic properties, 
mostly related to the recursions that they satisfy. In Section 4 we consider several pos- 
itive definite kernels on the A-dimensional unit ball of an infinite dimensional Hilbert 
space. In particular, we prove a basic property of the Szego type kernel studied in ||] 
by characterizing its Kolmogorov decomposition. In Section 5 we discuss the problem 
of recovering the representation from orthogonal polynomials and we prove a Favard 
type result. We plan a more detailed study of applications to multiscale systems in a 
sequel of this paper. 

l 



2. Preliminaries 



We briefly review several constructions of the tensor algebra and introduce the nec- 
essary notation. We also review the connection with displacement structure theory as 
established in 0. 

2.1. Tensor algebras. The tensor algebra over is defined by the algebraic direct 
sum 

r N = ® k > (c N r k , 

where (<C N )® k denotes the /c-fold tensor product of with itself. The addition is the 
componentwise addition and the multiplication is defined by juxtaposition: 

(x ®y) n = ^ x k® Vi- 

k+l=n 

If {e 1; . . . , e^} is the standard basis of C N , then {e^ (g) . . . ® | 1 < ii, . . . , % k < N} 
is an ortho normal basis of T/y. Let F^- be the unital free semigroup on N generators 
1, . . . ,N with lexicograhpic order -<. The empty word is the identity element and 
the length of the word a is denoted by \a\. The length of the empty word is 0. If 
a = %\ . . . if. then we write e a instead of ® ... ® , so that any element of 7/v can be 
uniquely written in the form x = '^2 ae¥ + c a e a) where only finitely many of the complex 
numbers c a are different from 0. 

Another construction of 7/v can be obtained as follows. Let S be a unital semigroup 
and denote by F (S) the set of functions : S — > C with the property that 0(s) ^ 
for only finitely many values of s. This set has a natural vector space structure and 
Bs — {$s I s G S} is a vector basis for F (S), where S s is the Kronecker symbol 
associated to s G S. Also F (S) is a unital associative algebra with respect to the 
product 

It is readily seen that -Fo(F^) is isomorphic to T^. Since each element in Fo(F^) can 
be uniquely written as a (finite) sum, = ^ CTeF + c a 5 a , the isomorphism is the linear 

extension $i of the mapping 5 a — > e a , o G F^. 

Another copy of the tensor algebra is given by the algebra V% of polynomials in 
A" noncommuting indeterminates X\, . . . , with complex coefficients. Each element 
P G V% can be uniquely written in the form P = J2cre¥+ c <jX a , with c a ^ for finitely 
many <r's and X a = X i± . . . X ik where a — i\ . . . i& G F^-. The linear extension $ 2 of 
the mapping 5 a — > X a , a G F^-, gives an isomorphism of T N with V%. 

Yet another copy of T/v inside the algebra of lower triangular operators allowed for 
the connection with displacement structure established in Thus, let S be a Hilbert 
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space and define: So 
(2.1) 



£ and for k > 1, 



£k — £, 



k-l 



£ 



k-l 



C k-V 



N terms 



For £ = C we have that can be identified with (C JV )' x,fc and Z/v is isomorphic to the 
algebra £^ of lower triangular operators T = [7y] £ £(©fc> Cfc) with the property 



(2-2) 



V 

iV terms 



for 2 < j, z, j > 1, and Tj = for all sufficiently large j's. The isomorphism is given 
by the map $3 defined as follows: let x = (xq, x\, . . . ) G 7/y (a; P G (C^)® 31 is the pth 
homogeneous component of 2); then x p = J2\a\= P c a e a anc ^ f° r 3 — 0> Tjo denotes the 
column matrix [co-jj^.. ., where "T" denotes the matrix transpose. Then T j0 = for 
all sufficiently large j's and we can define T G £(©fc> Cfc) by using ( |2.2j) . Finally, set 
$s(x) = T. 



£.1?. Displacement structure. We can now describe the displacement structure of the 
tensor algebra. We write this connection for C° N . Then it can be easily translated into 
any other realization of the tensor algebra. Let = [T^] G £(ffifc>oCfc) , k — 1, . . . , N, 
be isometries defined by the formulae: Tij = for i 7^ j '• + 1 and Tj+i^ is a block-column 
matrix consisting of N blocks of dimension dim Q, all of them zero except for the kth 
block which is the identity on C^. We have the following result noticed in 

Theorem 2.1. LetT G £% and define A = I - TT* . Then 



N 



(2.3) 

where 



A-Y,FkAF* = GJ n G\ 



k=l 



G 



1 T o 
T01 



and J 



11 



1 
-1 



The model C° N of the tensor algebra is also useful in order to extend this algebra to 
some topological tensor algebras - see [JT2"|. Here we consider only the norm topology 



and denote by Cn the algebra of all lower triangular operators T = [Tij] G £(©fc>oQ 
satisfying (|2.2|). 



2.3. Multiscale processes. Multiscale processes are stochastic processes indexed by 
nodes on a tree. They became quite popular lately, see fl|], f2|, and have potential to 
model the self-similarity of fractional Brownian motion leading to iterative algorithms 
in computer vision, remote sensing, etc. 

Here we restrict our attention to the case of the Cayley tree, in which each node has 
iV branches. The vertices of the Cayley tree are indexed by F^. 



Let (X, JF, P) be a probability space and let {v a } ae f+ C L 2 (P) be a family of random 
variables. Its covariance kernel is 



K(o, r) = v a v T dP, 
Jx 



and assume that the process is stationary in the sense (considered earlier, e.g. [|T0| 
that: 



(2.4) K(to, to') = K(o, a') , r, a, o' G F 



' r- H7+ 

JV' 



(2.5) K(o, t) = if there is no a G F^ such that o = ar or r = acr. 



Conversely, by the invariant Kolmogorov decomposition theorem, see e.g., [IS], Ch. II, 
there exists an isometric representation u of on a Hilbert space /C and a mapping 
f : F^r — > JC such that 

K(o,t) = (v(r),v(o)) K , 
u(t)v(o) = v(to), 

for all o, t G F^y-, the set {f(cr) | o G F^} is total in /C, and u(N) are 

isometries with orthogonal ranges. 

This class of multiscale processes would be suitable to model branching processes 
without "past". If a "past" should be attached to a process as above, we could try to 
consider processes indexed by the nodes of the tree associated to the free group on iV 
generators 1, . . . , N. As mentioned in the introduction, we plan to look at this matter 
in a sequel of this paper. Here we focus on processes with covariance kernel satisfying 
( p.4|) and ( |2.5|) . It was shown in || that such a kernel has displacement structure. 
Also, it is clear that for all j, k > 1, 

(2-6) [K(o,T)] ]aH!lTl=k = ([K(a , ,r J )] yH _ lM=k _ 1 )® N ) 

so that the kernel is determined by the elements s a = K(0, o), o G F^. 

By Theorem 1.5.3 in 0, each positive definite kernel K on F^ is uniquely determined 
by a family of contractions {7o-, r | o, r G F^, cr ^ r} such that j ajCT = 0, o G F^, and 
otherwise 7 CTjT G C(V ltT+1 T , D 7 * _ ) (for a contraction T between two Hilbert spaces 

D T = (I — y*T)V 2 denotes the defect operator of T and V T is the defect space of T 
defined as the closure of the range of Dt - note that in our case 7 CTjT are just complex 
numbers and the condition 7 CTiT G £(P 7ct+1>t , P 7 * ) for o -< r encodes the fact that 
|7<t+i,t| = 1 or |7(j |T _i| = 1 implies that 7^ = 0; also, r — 1 denotes the predecessor of 
r with respect to the lexicographic order -< on F^, while o + 1 denotes the successor 
of cr). In addition, the positive definite kernel K satisfies ( p.4|) and (|2.5|) if and only if 



(2-7) 7ra,ra' = 7<t,<t' , 7", C, (j' G F 



+ 



(2.8) 7(t,t = if there is no a G F^ such that o = ar or r = ever. 
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We define 7 CT = 70 j(J , a G F^ and we notice that {7^ \ a, r G F^, a ^ r} is uniquely 
determined by {7 -} (7gF + by the formula 



(2-9) M CTH ,| T | =& = (b^l^H-^i^.!)^, J, * > I- 

3. SZEGO POLYNOMIALS 

We introduce polynomials in several noncommuting variables orthogonal with respect 
to a positive definite kernel K satisfying Q2.4|) and (|2.5| ). We extend some elements of 
the Szego theory to this setting. 

The kernel K being given, we can introduce an inner product on .Fo(F^) in the usual 
manner: 

(3.1) M)k= K(a,r)i;(T)J{a). 

a,re¥+ 

By factoring out the subspace Mk = {0 G F (F^) | (</>,</>) # = 0} and completing 
with respect to the norm induced by ( |3.1| ) we obtain a Hilbert space denoted Hk- A 
similar structure can be introduced on V%. Let P = J2a&+ c aX a) Q = XLgf+ daX a 
be elements in V%, then define: 

(3.2) (P,Q) K = K (^ T ) d rCa. 

<t,t&+ 

By factoring out the subspace M.k = {P G V% \ (P,P)k = 0} and completing with 
respect to the norm induced by Q3.2|) we obtain a Hilbert space denoted L 2 {K). One 
can check that the map $ 2 defined by 5 a — > X a , a G F^, extends to a unitary operator 
from H K onto L 2 (K). 

From now on we assume that for any a G the matrix [K{cr, T)] aT<a is invertible. 
This implies that M.r — and V% can be viewed as a subspace of L 2 (K). Also, for any 
a G F^r, {X a } a ^ a is a linearly independent family in L 2 (K). Then, the Gram-Schmidt 
procedure gives a family {w} CTgF + of elements in V% such that 

(3.3) <f a = ^ aa, T X T , a a>a > 0; 

(3.4) {<p v ,(p T ) K = 0, 0^a^r. 

An explicit formula for the orthogonal polynomials ip a can be obtained in the same 
manner as in the classical (one variable) case. Define for a G F^-, 

(3.5) A, = det r')] CT ,, r , d(T 

and let {7<t} (j6 f+ be the parameters associated to K as described in Section 2.2.3. 
Note that since all the matrices [K(u,T)] aT ^ a , a G F^, are assumed to be invertible, 
it follows that | j a \ < 1 for all a G F^. 



Theorem 3.1. (1) (p$ = 1 and for -< a, 

1 



(3.6) 



fa 



det 



1 Xx ... X CT 



(2) For$^a = i!...i k , 
(3.7) <^ CT 



Yj^(X a + lower order terms). 



Proof. The proof is similar to the classical one. Thus, we deduce from the orthog- 
onality condition (|3.4j ) that {<p a ,X T i)K = for ^ t' -< or, which implies that 
J2 T ^a a v,TK{r' , t) = for ^ t' -< er. Using the Cramer rules for the system 

EhA^(t'.t) = 0, 0^r'^a, 
with unknowns a ff>T , we deduce 



<3-(7,(7 





■ [X(a', r') 




det 








1 Xi 





Therefore, 



Ax-i 



det 



I X l ... X a 

We now compute a CT)CT and -D^ in terms of the parameters {7o-} - e F+ of K. First we 
notice that 



(det 

1 Ii ... I, 
and since X a = -^—f a + c t^t, we deduce 



-Da = ( (/'a + > C t A t )k = — . 

n n < * n z 



^(7,(7 ^(7,(7 



T^<(7 



so that 



D„ 



^,(7 

which gives ( |3.6|) . 

In order to compute D a in terms of {ja}^^ we use Theorem 1.5.10 in || and the 
special structure of D a . Thus, 

d.= n (i-^vi 2 ) 
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and for -< a = i\ . . . we deduce 

1 D n 



n (i-K^r)- 

i<j<fc 



Then, 



IL 



<j<fc 



(l-|7i,...i fe P)V 



572(^(7 + lower order terms), 



which gives 



□ 



We illustrate this result for N = 2. From now on it is convenient to use the notation 

d CT =(l-| 7(T |2)V2 )(TGF +_ { } . 



Example. Let N = 2 and assume the positive kernel K satisfies the conditions in 
Theorem [Oj . We have = 1 and the next three determinants are: 

1 s 1 
si 1 



D 1 = det 



4; 



did: 



2^2. 

2) 



D u = det 



ill UoClii • 



1 Si s 2 
D 2 = det si 1 

_S 2 1 

1 Si s 2 s u 

Si 1 Si 
s 2 1 
In si 1 

Using Theorem |3J] we can easily calculate the first four orthogonal polynomials of K. 
Thus, ip% = 1 and then 

1 x x 



1 , 

— det 
di 



-— + —X 1 ; 



<^2 



d\d 2 



det 



JT2_ + 7i72 x 
dic?2 d\d 2 d 2 ' 



1 Si s 2 
si 1 
1 Xi X 2 

where we used the fact that s 2 = d\^ 2 . Then, after some calculations 

1 S\ S 2 Su 

Si 1 Si 
s 2 1 

1 X\ x 2 xf 



<P11 



d^d^dn 



det 



7n 



+ 



7i 



TiiTi s Y , TiiTa Y 1 ^2 



did 2 d\i d\di\ d\d 2 d\\ d 2 d\\ d\\d\ 

We establish that the orthogonal polynomials introduced above satisfy equations 
similar to the classical Szego difference equations. 



Theorem 3.2. The orthogonal polynomials satisfy the following recurrences: <p$ = 1 
and for k E {1, ... , N}, o E 



(3.8) 

where 

(3.9) 



"fccr 



1 and for k E {I,... ,N}, a E F^, 

1 



Proof. We deduce this result from similar formulae obtained for an arbitrary positive 
definite kernel. In this way we can show the meaning of the polynomials ip^ a , a E F^. 
Let [ti y j\ i j >l be a positive definite kernel on N and assume that each matrix A^'^ = 
[tk,i] i<k ;<„•, 1 < i < j, is invertible. Also, assume t k ,k = 1 for all k > 1. Let Fjj be the 
upper Cholesky factor of A™', so that i^^ is an upper triangular matrix with positive 
diagonal and A™' = F*jF it j. A dual, lower Cholesky factor is obtained as follows: 



define the symmetry of appropriate dimension, 

... 



J 





I 

1 



/ 

1 







and then let F it j denote the upper Cholesky factor of = JA^J. If G it j = 

JF id J, then 

A a,j) = ji r -rj = jFfjFijj = (■]/:, .; r 

and Gij is a lower triangular matrix with positive diagonal, called the lower Cholesky 
factor of A^\ Let Pij be the last column of F^ 1 and let P\- be the first column of 
Gjj, that is 



Pi,j = FrlE, P* j = G$JE, 

where E — [ ... / ] r . Let {^i,j}i<i<j be the parameters associated to [ti t j\ i ?>1 
by Theorem 1.5.3 in || and let pij 



[l - \r id \ 2 ) 1/2 . We have that 



(3.10) 
(3.11) 



Pl,n 
pi . 



1 







r l,n-\ 




P2,n 


dl,n 






"1,71 







+ 



«l,n 







These formulae are presumable known to the experts. For the sake of completeness 
we give a proof here based on results and notation from ||. First we introduce the 
following elements. For i < j, 



(3.12) Lf = L({r hk } 



\i=i+i) 



[ r i,i+l Pi,i+l r i,i+l 



Pi,i+1 ■ ■ ■ P 



1,3 



and 



C 



(0 



r i+1;j p i+ 2,j ■ ■ ■ Pj-ij 

. r i,jPi+l,j ■ ■ ■ Pj-lj 



K 



r i,i+lPi,i+2 ■ ■ ■ Pi,j 



fi,j-lPi,j 



Phj 



Also, we define inductively: D\ ' ' = Pi^+i, 



(3.13) 



D? = D({r i)k y k=i+1 ) 







'•j 



P 



We also need to review the factorization of unitary matrices. This is an extension of 
Euler's description of 50(3). First we define 



Rj~i{ r i,k) — Ik-1- 



T"i,k Pi,k 
Pi,k i,k 



where Ik-\~i is the identity matrix of size k — 1 — i. Then, 



R"i,j Rj—i^iji+l) ■ ■ ■ Rj—iij* 



1,3 ) i 



and 



It turns out that any unitary matrix can be written as a matrix of the form of Uij. 
The main idea for the proof of ( |3.10| ) is to use the identity 



(3.14) 



which follows from the relations (1.6.10), (6.3.8) and (6.3.9) in Thus, we notice 
that (|3TTJ) implies 



P* . = F U- E 

1 l,j 1 1,3^1,3-^1 



which is more tractable than the original definition of Pfj. This is seen from the 
following calculations. Using formula (1.5.7) in ||, the above definition of D^, and 



the notation D x ^ = (D± l> ) 1 , we obtain that 



Pl,n 



i -4 n) Dr (n) 

Lo *£i?r (B) 



1 Pl.n 1 







1 

Pl,n 



-(n-l) ri,„ n -(n-l) j^(n-l) 



Pl,n 



pi. 



(n) 



r 2,n 



Pi. 



1 

Pl,n 



Pl,n 1 1 
1 

Pl,n 



Pl,n 



Pl,n 





P2,n 



r%n jj-(n-l) Tv-(n-l) 

Pl.n 1 1 
1 



r 2,n 



Pl,n 1 



pi. 



(n-l)-l (n-l) 



+ 



r l, 



Pi,' 



• • Pl,n-1 



r 2,n 



D -(„-l) K (n-l) 



The proof of formula (1.6.15) in [Q gives 



t (n-l) n -(n-l) T^(n-l) . 



Pl,2 • • • Pl,n-1 



and using formula (1.5.6) in |J we deduce 



r 2,n 



^-(n-l)^-!) 



ci-l n -(n-l) T^(n-l) 
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therefore 



Pi, 





P2,n 



Pl,n 



1 

Pl,2—Pl,n-1 



p-l n -(n-l) Tfin-1) 
~ r 2,n-l- u l 



It remains to show that 



r l,n-l 



1 



Pl,2---Pl,n-1 



TTi-l n -(n-l) Tfin-i) 
~ r 2,n-\ u \ 



To that end we notice that using formula (1.5.8) in M, the definition of £/i )n -i, the fact 
that -Ri jn -i is a unitary matrix, and the notation £*( n_:L ) = we obtain 

i 







pi,2.--pi,n-l 

1 r *(n-l) 



^2,n-l Pl , 2 ...pi, n -i 2,n-l^l 



Pl,2 
1 



-^1,71-1-^1 



n-1 



C/ 2 ,n-l 
1 



PI, 2---P1, n-1 



P1.2..-P1, 



p-i rn 



It follows that all we have to show is the equality 



(n-1) jrin-l) 



rp—1 j * 



*(n— 1) 



p-i r>~ y ! K y 

P\,2 ■ ■ ■ P\,n-\ 

Now this is a simple consequence of the formula T*D?* = D?T* for the contraction 
T = Lf' 1] . The formula (|3TT| ) can be proved in a similar manner. 

We rewrite (|3.10|) and (|3.11|) for a positive definite kernel K satisfying (|2.4j) and 
( p.5|) . We notice that ff n _i is replaced by <^| (T _ 1 and then we have to show that P 2 ,n 
can be expressed in terms of <p a . This follows by taking into account the relations (|2.7|) 
and (|2.8| ) and using systematically ( 3.10|) . We can omit the details. □ 



The previous recurrence equations look quite similar to the classical Szego recursions, 
still they have a new component that is best understood when considering (|3.10|) and 
( p.ll|) . This type of recurrence equations was also found in 0] in connection with some 
derivations on Cn- It turns out that these derivations are related to those considered 



in 



11 , and later studied in 



We also notice a graded form of the recurences ( |3.8|) and (|3.9| ). It is convenient to 
introduce the following notation: for n > 1, we use ( |3.12| ) and define 



9n = L({7 <T }| CT | =n ). 

ll 



It was explained in || that g n , n > 1, are the parameters associated to the kernel K 
in ||lb| . We also use ( |3 . 1 3[ ) in order to introduce the notation 

H n = £>({7a}n=n), n > 1. 

Let cr(n) be the largest word (with respect to the lexicographic order) of lenght n, that 
is cr(n) = TV^JV. 

JV terms 

Corollary 3.3. TTie Szego polynomials satisfy the recurrences: for k > 1, 
(3.15) [</>]=fc] = ([ x i ■■■ x n] [p\a\=k-i]® N - ¥>i (fc )_i0*) *5 

( 3 - 16 ) ik 1 (- 1 ••• ^]h CT |= fc -i] eJv ^+^ (fc) -i)- 

|r|=fc 

Proof. Both ( P-15| ) and ( |3.16| ) follow by direct calculations from Theorem \!>.2\ □ 



4. Christoffel-Darboux formula 

A first consequence of the Szego formula in the classical case is the Christoffel- 
Darboux formula. Here we find a similar formula in several noncommuting variables. 
To that end we introduce additional notation. Let £ be a Hilbert space. In this paper 
£ will always be infinite dimensional. The iV-dimensional unit ball of £ is defined by 

B N (£) = {Z=(Z 1 ... Z N )\(Z\Z)<I S }, 

where for two elements Z = ( Z\ ... Z N ) and W = ( W\ . . . Wn ) in C(£) N 
we define 

N 

(4.1) {Z\W) = Y,Z k W* k . 

k=l 

We also need a sort of Szego kernel for B N (£). One suggestion was given in j|] to 
consider the following construction. For Z G B^(£) define 

(4.2) E{Z) = [Z a }™ ]=0 e C(® k > £ k , £). 

Also we use the notation diag(S') to denote the diagonal operator in C(@k>o£k) with 
diagonal S. A Szego type kernel on Bn{£) is given by the formula 

K S (Z,W) = E(Z)E(W)*, Z,WeB N (£). 

The next result explain two important properties of K. 

Lemma 4.1. (a) For any T e £{£) and Z,W G B N {£), 

N 

E{Z)diag{T - Z k TW;)E{W)* = T. 

k=l 

{b) The set {E(W)*£ | W G B N {£)} is total in ® k > £ k - 
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Proof, (a) Using directly the definitions, 

E(Z)di ag (T - Zti Z k TW* k )E{Wy 



T + E W >i Z a TW; - Ef =1 E(Z))dizg(Z k TWZ)E(W)* 



= T + E H >i W - ELi E k |> z„z k Tw;w; 

= T + E H >i W - E M >i Z a TW* 
= T. 

(b) Let e = {e -} cre F+ be an element of (B k >o8 k orthogonal to the linear span of 
{E(W)*8 | W E B N (8)}. Taking W = 0, we deduce that e = 0. Next, we claim that 
for each a G F 4 ^ — {0} there exist 

W l = {W[,... ,W l N ) eB N (8), 1 = 1,... ,2\a\, 

such that 



wi 1 ... w? wl 



range 
and 

W* = for all r^a, |r| > \a\, I = 1, . . . , 2\a\. 

Once this claim is proved, a simple inductive argument gives e = 0, so {E(W)*£ \ 
W G Bn{£)} is total in ® k >o8 k . Therefore we focus on the proof of the claim. 

Let {efj}fj = i be the matrix units of the algebra M n of n x n matrices. Each e™- is 
annxn matrix consisting of 1 in the (i,j)th entry and zeros elsewhere. For a Hilbert 
space £\ we define E^ = e™- <g> Ig 1 and we notice that 

(4.3) E-E-=6 Jk Ei /•;;• 

Let 8 be infinite dimensional, so that 8 = £® 2 ' CT ' for some Hilbert space 8\. Let 
a — i\ . . . i k . For s — 1, . . . ,N, we define 

J s = {I e {1, . . . , k} | = s} 

and 

W ^ = ^E £ vS-W s = l,...,iV, p=l,...,\a\. 

T^zJ s 

We show that for each p e {1, . . . , |cr|}, 



(4.5) W? = for r^a, |r| > \a\. 
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Using Q4.3|), we deduce 



2 Z_/s=l Z_/rgJ s r+p,r+p— l-'-V+p— l,r+p 
1 y^V 



n 2ff 



p2|<r| 
■<reJ s ^r+p,r+p 



hence iy p G B^(£) for each p 
deduce by induction that 



1 y^/c p2|o-| 

2 2—ir=\ r+p,r+p 

1, . . . , |<t|. For each word t = j\ ■ ■ ■ 3k £ 



{0} we 



(4.6) 



r+p+fc— 1> 



where A r = J ife _ p - p) C {1, . . . , N} and J jfc _ p - p = {I - p \ I G J ifc _J. 

We show that A a = {1} and A T = for r ^ a. Let g G A r . Therefore, for any 
p G {0, . . . , k — 1} we must have q + p G Jj fe or ik+i- q - p = For p = — 1 we 

deduce j\ = 12-q and since 2 — q > 1, it follows that q < 1. Also g > 1, therefore 
the only element that can be in v4 r is q = 1, in which case we must have r = a. 
Since Z G Ji k+1 _ { for each Z = 1, . . . , k — 1, we deduce that A CT = {1} and A T = for 
t ^ a. Formula ( |4.6|) implies (fOJ). In a similar manner we can construct a family 



P 



a + 1, 



2 1 a I, such that 



W, 



*p 



1 



2\a\ 
p+k,pi 



and 



Thus, for s — 1, 



W r p = for r ^ a, 
N, we define 



T > <7 . 



{Z G {1, . . . , k} I z fc = s} 



and 



7! ^ ^ 



r+p— fc,r+p— k— 1' 



p 



M + 1, 



2\a\ 



T&K S 



Now, 



w. 



*1 



*2|cr| 



2 fc L 



p 2|o-| 



p 2|<7| p 2|a| 



2kl 

2fc,fc 



whose range is £. This concludes the proof. 



□ 



We note that the result given by Lemma |4.1| (b) is not true in case £ is finite dimen- 
sional. The meaning of the result is that in case £ is infinite dimensional then E is 
really the Kolmogorov decomposition of the kernel K$- 
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We now let a polynomial P = £ CTeF + c o-X a G V% take values on Bn{£) by the 
formula 

(4.7) P{Z) = c ° z a, z e ®n{£). 

Define the Cristoffel-Darboux kernel by the formula 

(4.8) K CD (Z,W) = E(Z)dmg(ifl in) (Z)cpl in) (Wy - ^ ip T (Z)(p T (W)*)E(W)*, 

\r\=n 

for Z,W G B N {£). 

Theorem 4.2. For any Z,W E B N {£), 

k cd (z,w)= MZ)Mwy. 

0<\r\<n 

Proof. From (|3.8|) and (|3.9|) we deduce 

vLWvLW - <Pk*Wv>ka(wy = vL-iWvL-i(wy - z k¥ > ka (z) Vka (wyw^ 

for any k G {1, . . . , A^}, a G F^, and Z,W G Bn(£). Adding all these relations for 
k G {1, . . . , N} and < \a\ < n — 1, we deduce 

iV 

0<\a\<n fc=l 0<|cr|<n-l 

give 



This relation and Lemma [O 
K CD (Z,W) 



E(z)(Eo<w\< n Mz)Mwy - ELi ^(£ < H <„ WW) W 
E Q <w<„^(^)^W*. 



□ 



We can show one more application of Lemma [4.1| . For a formal power series 

/ = c<jX<J ' 

<x G F+ 

in A" noncommuting variables ATi, . . . , X^r, we denote by Tf the lower triangular infinite 
matrix associated to / as described in Section 2.1. We denote by Sn the Schur class of 
those formal power series / with the property that Tf is a contraction in £(©fc> Cfc). 
If £ is an infinite dimensional Hilbert space then we can define f(Z) for Z G Bn{£) as 
in by the formula 

(4.9) f(Z) = E(Z)(T f ®I £ )/£. 

We notice that this definition is consistent with ([4.7] ) . We extend a familiar character- 
ization of the Schur class to the setting of this paper. 
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Theorem 4.3. The formal power series f belongs to Sn if and only if 

C f (Z, W) = E{Z)diag{I - f{Z)f{Wf)E{W)\ Z,W e B N (S), 
is a positive definite kernel on B^{£). 

Proof. Using Lemma 3.1 in we deduce that for Z,W G Bjv(£ ), 
E(Z)(I - (7) <g> I £ )(T f ® I e )*)E(W)* 

= E(Z)E(W)* - E(Z)(T f ® I £ )(T f ® J £ )*)£(H0* 

= - £(Z)diag(/(Z))diag(/(WO*)£W 

= £(Z)diag(7 - = C/(Z, W). 

This relation implies that if / G <S/v then C/ is a positive definite kernel on Bn(S). For 
the converse implication we have to use in addition Lemma |4. 1|. □ 



5. Inverse problems 

In this brief section we prove a Favard type result for orthogonal polynomials in 
several noncommuting variables. 



Theorem 5.1. Let {7o-} (je F+ be a family of complex numbers with 70 = and \ < 1 
for a G F^r — {0}. Then there exists a unique positive definite kernel K satisfying (|2.4j) 
and (|2.5| ) such that the polynomials <p a , a G defined by the recursions: <p$ = 1, 



1, and for k G {1, . . . , N}, a G F^, 

1 

d k 



Vka = —(X k ip a - Tfc^L-l •■ 

a k 

Uka 

are orthogonal with respect to K . 

Proof. Once again we rely on some results that are known for arbitrary positive definite 
kernels on the set of integers. In this way, the proof is quite straightforward. Let 
{7 CTjT I a,r G F^r, a z< r} be the family of complex numbers associated to {7 -} - eI p+ by 

( p.9|) . Let K be the positive definite kernel associated to {-f a ,r \ o~, r G F^,a ^ r} by 
Theorem 1.5.3 in j3|. By Theorem |3~2| , the polynomials <p a , a G F^, defined by the 
recurrences: (p$ = 1, ipt = 1, and for k G {1, . . . , N}, a G F^, 



1 



Vka = -j—(X k ip a ~ ^kaf i 



ka 



a-\h 



<PL = -l-(-lka X k<fa + V>L-l)> 
Uka 

must be the orthogonal polynomials of K. □ 
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